Study of the Hadronic Transitions T(2S)^ (77, tt°)T(1S) at Belle 
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We study the rare hadronic transitions T(2S) — > T(lS)r; and T(2S) — ¥ T(lS)n° using a sample 
of 158 xlO 6 T(2S) decays collected with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. We measure the ratios of branching fractions (3) • B( ^^g^x(is]n+'l- ) ~ (1.99±0.14 (stat) 

±0.11 (syst)) xlO -3 and B fr[2S)ir(is)J+l- ) < 2 ' 3 x 10-4 at the 90% confidence level ( CL )- 
Assuming the value S(T(2S) ->■ T (1S>-7T + ) = (17.92±0.26)%, we obtain S(T(25') -> T(lS)r?) = 
(3.57 ± 0.25(stat) ±0.21(syst)) x 10 -4 and B(T(25) T(1SV°) < 4.1 x 10" 5 (90% CL). 

PACS numbers: 14.40.Pq,13.25.Gv 



In recent years, hadronic transitions between quarko- 
nia have led to an impressive series of discoveries pQ: 
AT(3872), Y"(4260), as well as h c and hi, were observed 
in transitions either from or to the ip and T states. 
The phenomenology of these transitions is commonly de- 
scribed with the QCD Multipole Expansion formalism 
(QCDME) which allows one to classify the tran- 

sitions in a series of chromoelectric and chromomagnetic 
multiplets. In particular, theoretical predictions for n 
and 7T° transitions [H |S] among states are being chal- 
lenged by experimental measurements. The r\ and ir° 
transitions between vector bottomonia should be medi- 
ated either by two Ml gluons or by one El and one M2 
gluon: both cases imply a spin flip of the b quark. The 
corresponding amplitude should scale as 1 /nib, an d its 
measurement yields information about the chromomag- 
netic moment of the b quark. 

By scaling from the ip{1S) — > J/ipn transition, one 
expects a transition width of: r[T(2S f ) — > T(lS)n] = 
0.0025 x T[ip(2S) J/ijj'n], and therefore a ratio of 
branching fractions 3< v<w + w - = 3( T (2S)->T(is)7r+7r-) ~ 
4.4 x 10 -3 [6]. From |4], one can calculate 3?,,^+^- = 
2.3 x 10~ 3 assuming the b quark mass to be = 
4.67 GeV/c 2 [7]; the prediction from [5] is Ji rit7r + 7r - = 
1.7 x 10~ 3 . A further suppression is expected for the ir° 
transition, which violates isospin; here, theory predicts 
r[T(25) ^ T(lS*)7r°] = 0.16 x T[T(25) ->■ T(15)r?] . 
The T(2,3,45') — > T(lS)ri transitions have been studied 
by BaBar 8 ( 9 and CLEO [B]; the measured branching 



fractions are either unexpectedly large (T(45)) or too 
small (T(2S*) and T(3S')). The parameters of the quark 
wave functions must be changed by more than 15% in or- 
der to account for these discrepancies [TU]- Searches for 
the 7T° transitions have only yielded upper limits [BJ [S] . 

We report here a new measurement of the transition 
T(25) ->■ T(lS)ry and a search for T(2S) ^ T(15*)7r° 
using the Belle detector at the KEKB c + c~ collider |11) . 
The T(IS') is reconstructed in both the e + e~ and ^ fi~ 
decay modes; we reconstruct the r\ meson via its decay to 
77 and 7r _ 7r + 7r°, and the ir° in the 77 final state. As a 
normalization sample, we reconstruct the dominant tran- 
sition T(2S') — > T(1S)tt~tt + , which has a branching frac- 
tion of (17.92±0.26)% [7j. The data sample for this anal- 
ysis includes an integrated luminosity of 24.7 fb _1 at the 
T(2S') resonance peak, corresponding to (158±4)xl0 6 
T(2S') decays, and an additional 1.7 fb _1 at y / s=9.993 
GeV used to study the QED continuum backgrounds. 

The Belle detector is described in detail elsewhere [H] . 
Here, we summarize the features that are relevant to the 
current analysis. The momentum of each charged track is 
measured using a four-layer double-sided Silicon Vertex 
Detector and a 50 layer Central Drift Chamber, embed- 
ded in a 1.5 T magnetic field, which allows tracking of 
charged particles with transverse momentum as low as 
50 MeV/c. Electrons and photons are detected in a large 
array of CsI(Tl) crystals (ECL) also located inside the 
solenoid coil. An iron flux return located outside the 
coil is instrumented to detect mesons and to identify 
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muons. 

Monte Carlo (MC) samples of the signal and of the 
dominant peaking backgrounds were generated using 
EvtGen |13j . Dipion transitions were generated assuming 
that the amplitude is dominated by the S'-wave contri- 
bution. QED continuum processes e + e~ — > (e + e~ or 
/j + /j~) + nj were generated using KKMC [T3]. The de- 
tector response was simulated using GEANT3 [T5], and 
beam backgrounds were accounted for using random trig- 
gers taken during each period of data taking. Final state 
radiation effects are accounted for by using PHOTOS [16] 
in EvtGen simulations. 

Charged tracks with momentum p* in the center-of- 
mass frame of the colliding e + e~ pair (CM frame) greater 
than 4 GeV/c are selected as candidate leptons from 
T(15) decay. In the following text, all the quantities 
computed in the CM frame are denoted with an asterisk. 
Electrons and muons are identified by the ratio i? e , M be- 
tween the electron (muon) hypothesis likelihood and the 
hadronic hypothesis likelihood. A track is identified as a 
lepton if R e or i? M is above a threshold value of 0.2 and 
then as a muon if > R e . The identification efficiency 
is 93.2% for electrons and 92.6% for muons. Pairs of 
lepton candidates with opposite charge and an invariant 
mass in the range 9.0 GeV/c 2 < M(£+£') < 9.8 GeV/c 2 
are then selected for further analysis. 

In order to reduce the effect of final state radiation 
(FSR) and bremsstrahlung, the momentum of all photons 
detected in the ECL within 200 mrad of each leptonic 
track is added to its momentum. A mass-constrained 
kinematic fit performed on the T(IS') candidate lepton 
pair is required to have a confidence level CL^s > 10~ 5 . 

A requirement on the polar angle in the CM frame 
of the e~ track with respect to the beam direction, 
cos(6**_) < 0.5, is imposed in the reconstruction of the 
V 77»T(LS) — * e + e~ final state in order to sup- 
press singly- or doubly-radiative Bhabha events, which 
represent the dominant QED background for this chan- 
nel. The Bhabha requirement is not included in the 
T(2S) —> T(1S)tt° analysis since the CLig requirement 
provides sufficient suppression. 

Dipion candidates used both for ij — > ir + ir~ir° and the 
T(2S) — > T(lS')7r + 7r _ transition are composed of oppo- 
sitely charged pairs of tracks, each of which has a distance 
of closest approach of less than 1.5 cm (0.5 cm) in the ax- 
ial (transverse) direction relative to the beam line. The 
cosine of the angle in the CM frame between these tracks 
is required to be less than 0.6 in order to reject events 
with photons that convert in the inner detectors. 

In the 77 — > 7T + 7r~7r° analysis, each photon produced 
in the ir° decay must have E\ a b > 57 MeV, the optimal 
threshold for the rejection of photons arising from beam 
background, and E* < 220 MeV. The n° candidate is 
then selected as the 77 pair with invariant mass closest 
to the nominal ir° mass [7]. The threshold values are 
optimized by maximizing the figure of merit, defined as 



FoM — -7=pjj where s is the signal yield and b is the 
background contribution. 

In the 77 — > 77 analysis, photons with 180 MeV < 
E* < 360 MeV are subjected to the requirement on the 
opening angle in the CM frame of cos0* 7 < —0.88 in 
order to reject combinatorial background. Events with 
more than one 77 candidate that satisfies these conditions 
are found to be a negligible fraction of the total Monte 
Carlo sample, and are rejected without introducing any 
further selection. 

The planarity of the event is exploited in order to se- 
lect the T(25) -> T(lS)iT decay. We select the pair 
of photons with CM momentum p* 7 that minimizes the 
scalar product (p* 7 • U£+£-)/|p* 7 | (where Ug+g- is a vec- 
tor normal to the plane formed by the dilepton pair in 
the CM frame) as the tt — > 77 candidate. 

After 77 or tt° selection, the T(liS') and 77 or tt + tt^tt° 
are subject to a second kinematic fit, constraining them 
to have the T(2S) invariant mass. The minimum confi- 
dence level CL2S is optimized for each decay mode using 
MC samples. For the ?y — > 77, CL2S is required to be 
greater than 6 x 10~ 4 (2 x 10 -3 ) for fi + (e + e~) events; 
for -q — > 7T + 7r~7r , the event is accepted if the fit con- 
verges. Finally, in the analysis of T(25) — > T(1S)tt°, 
CL25 must be greater than 10 -5 for both T(IS') decay 
modes. 

Since the signal events are fully reconstructed and the 
total momentum of all the charged tracks and photons 
in the CM frame is expected to be close to zero, a se- 
lection on p* tot = |pt(is) +f>7j,7r | is imposed, requiring 
p* tot <0.07 GeV/c in T(15) A e+e~,r) -> 77 events, 
which are more contaminated by radiative Bhabha, and 
Ptoi < ^-l GeV/c when investigating the other modes. 

The requirements described above result in an almost 
complete rejection of the QED and T(25) -> T(lS')7r 7r 
backgrounds. 

The T(2S') -> XbjQPhi -> T(lS')7i7 2 decay has the 
same event topology as 77 — > 77 and 7r° transitions. The 
kinematic limit of the less energetic photon, which arises 
from T(2S) Xbj0-P)li decay, is 162 MeV in the CM 
frame. Therefore, this background is completely rejected 
for 77 —¥ 77 by the photon energy requirement mentioned 
earlier. In the search for the ir° transition, this back- 
ground is still larger than expected signal; hence, we re- 
quire the energy of the less energetic photon to be higher 
than 170 MeV. This requirement rejects 99.5% of the 
XbjO-P) background and retains 34% of the signal events. 

The 7r + 7r~ transition represents a significant source of 
background only for the r\ — > 7r + 7r~7r° channel; in this 
case, we require AM = M(tt + tt~£ + £') - M(£ + £~) < 
0.44GeV/c 2 , since this observable peaks at AM = 
M(T(25)) - M(T(1S)) = 0.56 GeV/c 2 in tt+tt- events; 
this cut rejects 99.92% of the background and retains 
99.1% of the signal. 

Signal efficiencies for the various final states are sum- 
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marized in Table ||] The efficiency in the e + e~77 mode 
significantly differs from the one in /z + /x _ 77 since this 
channel is affected by a Bhabha veto included at the trig- 
ger level. A trigger simulation is used in order to account 
for this effect. 

The signal yield is extracted with a simultaneous, un- 
binned likelihood fit of the r) mass distribution in four 
different final states with a common branching fraction, 
as shown in Fig. [T] For the r\ invariant mass peak, we 
use a double Gaussian with parameters, that differ from 
channel to channel and fixed at values determined by the 
simulation. The background probability density function 
(PDF) shape, a Crystal Ball function [17] in r\ — > 77 and 
a Gaussian shape in r\ — > 7r + 7r~7r°, is chosen using the 
MC simulation; the parameters of the chosen PDFs, in- 
cluding the background yields in each channel, are left 
free in the fit of the four final states. The sum of the in- 
variant mass distributions for the four independent final 
states is shown in Fig. [2] 



The ratio of the branching fraction for the 77 transition 
to that for the dipion transition is given by: 



B(T(15)»7) _ N nJ 



At 



r),7T' r 7T'~ 



S(T(lS)7r+7r-) AT« Sfa-*/) 



e v,f 



where the r\ final state / is 77 or ir + ir ir° , N v is the sig- 
nal yield from the fit and is the number of detected 
T(2S) -> T(15)7r+7r- transitions for each T{1S) ->■ l+l~ 
decay mode: = 228167, N%% = 276261 and = 
504428. The charged dipion transition is selected with 
the same cuts as in the r\ — > 7r + 7r~7r° selection and requir- 
ing the event to have exactly one dipion and no ir° or 77 
candidates. According to the Montecarlo simulation the 
background contribution is negligible, and the efficien- 
cies for the normalization channel are e^. =31.26% and 
e££ =37.94%. The number of T(2S*) T(1S> events 



extracted from the simultaneous fit is N„ 



241 ± 17. 



The resulting ratio is = (1.99 ± 0.14) x 10~ 3 , 

where the error is statistical. 

In order to study systematic effects, the ratios 0i n 
measured separately in four different final states are re- 
ported in Table [n] The fitting procedure used for the 
individual channels is the same as that used for the fit to 
the full sample. 

In the T(25) —> T(lS')7r analysis, the 77 invariant 
mass distribution is fitted to a Gaussian function for the 
signal and a third-order polynomial for the background 
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TABLE I: Signal efficiencies 


15 






T(15) -»■ e+e" T(1S) -> 
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> 77 8.4% 25.7% 




n - 


► -K+n-TT 6.4% 7.6% 
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7T° 


-> 77 6.0% 7.8% 
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FIG. 1: Mass distribution of the 77 candidates in 
(top), 77e + e~ (middle), -k + n~ -k° fj, + n~ , and tv + tt~ 7r° e + e~ 
(bottom) final states. The fit function in blue represents the 
simultaneous fit of the four channels; the red curve shows the 
best-fit background component 
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TABLE II: Ratio % v 
pies. 



extracted from different subsam- 



Final state 



77^ 2 7 T(1S) 

77^ 2 7 T(15) - 

V ^3tt T(15) 

77^3tt T(15) 



2.16 ±0.19 (stat.) 
e + e" 2.15 ±0.38 (stat.) 
/x+jt/T 1.66 ±0.39 (stat.) 
e + e" 1.31 ±0.56 (stat.) 



Simultaneous fit 



1.99 ±0.14 (stat.) 




FIG. 2: 77/7r + 7r 7T° invariant mass for T(2S) 
candidates, summing all the four final states. 
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FIG. 3: Final fit to the 77 invariant mass for T(2S) 
T(1S)7T° candidates. 



TABLE III: Sources of systematic uncertainties. 



Source 


77 channel 


7r° channel 


Bhabha veto 


±2.5% 




Kinematic Fit 


±1.5% 


±1.5% 


Background fit 


±2.1% 


±6.3% 


77/71-° reconstruction 


±4% 


±4% 


Signal PDF 


±2% 


±1.8% 


Total 


±5.7% 


±7.8% 



(Fig. [3]). All parameters are freely varied except the 
width and the mean of the Gaussian, which are set to the 
values determined by MC simulation. No clear evidence 
for a 7T° signal is found in either the T(1S) — > e + e~ or 
the T(IS') — > mode. The signal yield from the fit 

is N„a = 10 ± 5. 

An upper limit on the number of T(lS)ir° candi- 
dates, N^ , ^ s determined by generating 5000 pseudo- 
experiments for different values of the signal yield rang- 
ing from to 30, using a Bayesian-frequentist hybrid ap- 
proach [TS] based on the ratio of CLs between the p- 
value of the signal±background hypothesis and the p- 
value of the background-only hypothesis. The resulting 
upper limit is N^ =20.5. 

ihe upper limit on the ratio s(T(2S)4-T(is)7r+7r-) 1S 
then calculated from the relation: 

S(T(1S)tt°) N^o 
*°>«+*- V(T(lS)n+n-) NH H(n° -+ 77) ■ e^ n 

where e^ is the efficiency reported in Table |lj reduced 
by a factor l-<7 sys =0.922 to account for systematic uncer- 
tainties. The contributions from each source, reported in 



Table |III[ are summed in quadrature in order to obtain 
the final uncertainty in each channel. The resulting up- 
per limit is 3? w o „.+„.- < 2.3 x 10~ 4 at 90% confidence 
level. 

The systematic uncertainty arising from the Bhabha 
veto procedure is obtained comparing by the signal yields 
obtained with and without the veto. The systematic un- 
certainty introduced by the 2S kinematic fit is studied us- 
ing a sample of T(25) — > Xbi,2l — > T(1S')77 events that 
are identified by modifying the r\ — > 77 selection. To re- 
construct these events, one photon in the range 90 MeV 
< E* < 180 MeV and one with 380 MeV < E* < 700 
MeV are required. Three different threshold values for 
CL25 have been used to estimate the systematic uncer- 
tainties on S(T(2S) X6i,27 T(15) 77 ). The sys- 
tematic uncertainty due to the choice of the parameters 
describing the background PDF is estimated by varying 
each within ±ler from the MC value in the 77 transition 
and changing the order of the polynomial fit in the ir° 
channel. The uncertainties arising from the possible dif- 
ference of signal PDF parameters between data and MC 
simulation are estimated by varying them within the er- 
rors and then comparing the obtained branching frac- 
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tions. The uncertainties related to the track reconstruc- 
tion and the total luminosity are canceled by the normal- 
ization to the T(2S) — > T(lS)ir + ir~ transition. An ad- 
ditional relative uncertainty of 4% due to neutral meson 
reconstruction is included. This error is determined from 
the discrepancy between data and MC in D° — » K~ir + Tr Q 
decay. 

In summary, using 24.7 fb _1 of data taken at the 
T(25) resonance peak energy, a measurement of the ratio 
3l„ j7r + jr - is obtained: 



(1.99 ± 0.14(stat) ± O.ll(syst)) x 10 
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This result is about 17% greater than the prediction of 
Ref . [5], 14% below the value extracted from Ref. [4] 
and less than half the value predicted by scaling from 
the ip{2S) — ¥ J/ipi] branching fraction [BJ. Assuming the 
branching fraction B(T(2S) T (lS>+7r-) = (17.92 ± 
0.26)% [7j, a new measurement of B(T(2S) ?(lS)rf) 
is obtained: 

S(T(25) -> T(lS)r)) = (3.57±0.25(stat) ±0.21(syst))xl0" 

where an additional systematic error of 1.4% is in- 
troduced in order to account for the uncertainties on 
S(T(25) ->■ T(lS')7r+7r-). This result is higher by about 
two standard deviations and more precise than those ob- 
tained by BaBar [9], and CLEO [6j. In addition, an up- 
per limit for the 52,0 ,+,- ratio, a factor of four more 
stringent than that of CLEO [BJ, is obtained: 

^o i7r+7r - < 2.3 x 10~ 4 , 

corresponding to the upper limit of 

S(T(25) -> T(15)7T°) < 4.1 x 10" 5 (90% C.L.). 

The upper limit for the ratio 

2(T(2S) -> T(15)tt°) 



S(T(25) -> T(lS)r)) 



< 0.13 (90% C.L.) 
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